Ultra-cold fermions in real or fictitious magnetic fields: 
The BCS-BEC evolution and the type-I— type-II transition 
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We study ultra-cold neutral fermion superfluids in the presence of fictitious magnetic fields, as well 
as charged fermion superfluids in the presence of real magnetic fields. Charged fermion superfluids 
undergo a phase transition from type-I to type-II superfluidity, where the magnetic properties of 
the superfluid change from being a perfect diamagnet without vortices to a partial diamagnet with 
the emergence of the Abrikosov vortex lattice. The transition from type-I to type-II superfluidity 
is tunned by changing the scattering parameter (interaction) for fixed density. We also find that 
neutral fermion superfiuids such as ®Li and ^''K are extreme type-II superfiuids, and that they 
are more robust to the penetration of a fictitious magnetic field in the BCS-BEC crossover region 
near unitarity, where the critical fictitious magnetic field reaches a maximum as a function of the 
scattering parameter (interaction). 
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A key experiment in the verification that neutral 
Fermi superfluids can evolve from the Bardeen-Cooper- 
SchriefFer (BCS) to the Bose-Einstein condensation 
(BEC) regime was the observation of quantized vortices 
throughout the BCS-BEC evolution [J upon rotation of 
the atomic cloud. This observation had a very dramatic 
impact beyond the atomic physics community, because it 
showed that superfluidity of Cooper pairs and of tightly 
bound bosonic molecules for s-wave pairing are the man- 
ifestation of the same type of physics. The key tool that 
permitted such realization is the tunability of the inter- 
action between fermions through the use of Feshbach res- 
onances. The same kind of tunability does not exist in 
■^He, the standard condensed matter neutral superfluid, 
or in superconductors. The situation is even worse in 
neutron and proton superfluids, which are thought to ex- 
ist in the core of neutron stars. 

Very recently, a new technique was developed that per- 
mitted the production of fictitious magnetic fields which 
can couple to neutral bosonic atoms [1, [1] . These ficti- 
tious magnetic fields are produced through an all optical 
Raman process, couple to a fictitious charge, but produce 
real effects like the creation of vortices in the superfluid 
state of bosons. In principle, the same technique can be 
applied to ultra-cold fermions, which coupled with the 
control over interaction using Feshbach resonances allows 
the exploration of superfluidity not only as a function of 
interaction, but also as a function of fictitious magnetic 
field. It is in anticipation of similar experiments involving 
ultracold fermions that we address in this manuscript the 
effects of fictitious magnetic fields on fermion superfiuids 
as a function of interaction. 

Unlike neutral superfluids, standard condensed mat- 
ter charged superfiuids (superconductors) can be of two 
types Many superconductors are now known to be 
type-II (including heavy fermions, organics, and high- Tc 
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FIG. 1. Universal phase diagram of the dimensionless fermion 
density nr^ versus scattering parameter l/{kFas), where = 
(f /{mc^) is the classical radius of a fermion with mass m and 
charge q\ cq is the speed of light, fc^ is the Fermi momentum 
and Os is the scattering length. The dotted line separates 
regions of type-I and type-II superfluidity. 



cuprates) , where the application of an external magnetic 
field beyond the lower critical field Hc^ leads to a non- 
uniform supcrfiuid phase, which appears in the form of 
the Abrikosov vortex lattice, until a second critical field 
is reached, when the system becomes normal. Other 
charged superfluids are known to be type-I and do not 
allow the magnetic field to penetrate the sample. These 
systems are perfect diamagnets until the critical field He 
is reached, where the charged superfluid becomes normal. 
The parameter that characterizes the type of charged su- 
perfluid is the Ginzburg-Landau parameter k = A/^ cor- 
responding to the ratio between the penetration depth 
A of the magnetic flcld into the sample and the coher- 
ence length ^ of the charged superfluid, such that type-I 
superfluids have n < l/-\/2 and type-II have k > 1/V2. 

In this manuscript, we study neutral fermion super- 
fiuids in the presence of fictitious magnetic fields and 
charged fermion superfiuids in the presence of real mag- 
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netic fields function of interaction (scattering pa- 
rameter). We show that throughout the crossover region 
between BCS and BEC superfluidity both '^Li and 
arc extreme type-II superfluids, and for charged super- 
fluids we find a phase transition from type-I to type-II 
superfluidity for fermions of density n = fc|./(37r^) inter- 
acting via a contact potential characterized by the inter- 
action parameter l/(fci?as). As shown in Fig.[Tl the phase 
boundary in the density n versus interaction parameter 
l/{kpas) occurs when the critical value Kc = 1/V2 is 
crossed. In the literature of charged superfluids the tran- 
sition from type-I to type-II was thought possible when 
induced by disorder and was described microscopically 
only in the BCS limit 0-13 • In contrast, here we show 
that, microscopically, a clean (no disorder) charged su- 
perfluid can exhibit a type-I to type-II transition induced 
by interactions. The phase diagram shown in Fig. [T]has 
a wider applicability to include standard charged super- 
fluids (like superconductors of condensed matter physics) 
and even proton superfluidity in nuclei or neutral stars, 
as long as the interactions can be described by a con- 
tact potential with corresponding scattering length a^. 
In addition, we indicate that neutral (charged) superflu- 
ids are more robust to the penetration of ficitious (real) 
magnetic fields near unitarity, where the critical fictitious 
(real) magnetic fields reach a maximum as a function of 
the scattering parameter. 

To describe the transition from type-I to type-II super- 
fluidity as a function of the interaction parameter and the 
properties of neutral (charged) superfluids in the presence 
of fictitious (real) magnetic fields during the BCS-BEC 
evolution for s-wave superfiuids in three dimensions, we 
start with the Hamiltonian density 



at the critical temperature Tc, where the order parameter 
vanishes. Here is the sample volume. The interaction 
g can be written in terms of the scattering length lead- 
ing to 1/g = -ml {\T,r?as) + (l/i')Ek [l/(2ek)]- The 
second self-consistency relation is the number equation 
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where f(^k) is the Fermi function, and 



g 



1 1 1 - f(ek+q/2) - f(e-k+q/2) 



^ Ck+q/2 



^-k+q/2 



ihuj 



(5) 



is the pair propagator, and uj is the Matsubara fi'cqucncy 
for bosons. 

The effective action is TS^s/fi ^ Eg ^"H9)|A(g)P + 
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A(9i)A*(92)A(g3)A*(gi-92 + 93) in terms 
of the order parameter A(g), where q = (q, io;). To 
study thermodynamic properties, we take itiuj = 0, or 
equivalently, A(r, r) = A(r), leading to the effective La- 
grangian density 
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Using the notation = tanh[^k/(2r)] and Ik = 
sech^[^k/(2r)], the coefficients of the Lagrangian density 
are a{^, ^ jj ^ Ek ^ ^^'^ constant term, 
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H(r) = Ev4(r) 
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where U{r) = / dr'y(r, r')V't(i")'/'j(r')^;(i"')'/'T(r) con- 
tains the attractive contact interaction potential 
V{r, r') = —g5{r~r'), and is the creation operator 

of fermions with mass m and spin a. Notice that g has di- 
mensions of energy times volume. To make progress, we 
rewrite the Hamiltonian H = J dr'H{r) from real space 
to momentum space 



■H = E ^kV'k.aV'k.a - 5 E ^k,q^k',q, 



(2) 
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where ^ 



k+q/2,t^-k-|-q/2,4. 



creates a fermion pair 
with center of mass momentum q and relative momentum 
2k, ^k = Ck — M is the kinetic energy term with ek = 
fi^k'^ /{2m) and yu is the chemical potential. 

Integration over the fermion fields [qI leads to the order 
parameter equation 



- = -y 



tanh [a/(2Tc)] 
2a 



for the coefficient of the gradient terms ViAV^A, and 
L^b{nc,Tc) = Ek - ^^'^ coefficient of the 

non-linear quartic term. Notice that Cij = c5ij for s-wave 
superfluids. 

In general, near Tc, a{fi,T) = — aoe(r), where 
ao = Tc [da/dT]r^^ and e{T) = (1 - T/T^). In the 
BCS limit of l/{kFCLs) —oo, L^ao = Pf, where 

the density of single parti- 
cle states per spin channel at the Fermi energy ep- 
Also the coefficient of the quartic term is Li^h = 
[7C(3)/(87r2r2)] Vp, while the coefficient of the gradi- 
ent term is L^c = [l (,{?,) /{12tt'^T'^)\V pep. Here, the 
zeta function C(3) = 1.202, while the critical temperature 
Tc = (8e''-V7r) ep e^]i[-TT / {2kp\as\)] , with « 1.781, 
and the chemical potential ^ — ep. However, in the 
BEC limit of l/{kpas) +oo, L'^a^ = Vptp / 
Correspondingly the coefficient of the quartic term is 
V^b = (7r/32) T> p / ep\ii\) , and the coefficient of the 
gradient term is L^c = [tt / Id)!) p / ep\ii\. In this case. 



(3) Tc ~ 0.218eF and ^l = Eb/2, where Eb = /{maj 
the two-particle binding energy in vacuum. 
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Next, we scale the order parameter to = ■y/cA(r) 
and introduce an external (real or fictitious) magnetic 
field via the vector potential A(r), using the substitution 
Vi — >■ Vi — 2iqAi/{hco), where q is the real or fictitious 
particle charge and cq is the speed of light. The difference 
in free energy density between the charged superfluid and 
its normal state in the presence of magnetic fields takes 
the Ginzbug-Landau form 
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where H = V x A is the real or fictitious microscopic 
magnetic field. The parameter a ~ a/c changes sign at 
T = Tc, however /? = b/c^ is always positive guaranteeing 
the stability of the theory. It is also useful to define the 
flux quantum (f>q = Trhco/q, which will be used below. 

Minimization of J-q l with respect to V-' and A lead the 
order parameter equation 



2m 



-V-I^A 
hco 
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and to the current density 



j = -^(V*VV^-^VV'*)-^|^pA. (7) 
im mcQ 

Using the relation V x H = 47rj/co and taking the curl of 
the current density leads to the London equation A^V x 
(V X H) + H = where A = ^?ncg/(167rq2|^|2) ig the 
magnetic penetration depth. Since \ip\'^ = |q;|//3 = |a|c/6 
in weak magnetic fields, the penetration depth becomes 
A(r) = AGL|e(r)|-i/2, where Agl = VV(167rr,aoc). 
Here, rq = q^ /(mc^) is the classical radius of a fcrmion 
with mass m and charge q in CGS units. Since jV'P 
plays the role of the superfluid density n^, we may 



write 



\MMT)\ = 



sfi\e{T)\, where 



"^sfl — IV'oP ~ ciQc/b is the temperature independent 
pref actor. This observation allows us to write kp\Qi^ = 
yJliiT / {l^kprq)] {n/risfi). The prefactor reflects a 
zero temperature extrapolation of the superfluid density 
rig, however, in a Galilean invariant system we must have 
ns,o ~ n/2. Indeed, in the BCS limit ris,o = n/2 such 
that kpXcL = ■\/37r/ (Skprq), while in the BEG limit 
rig^o = 3n/8 leads to kpXoL = \/ 1^ j i^kpTq). The BCS 
value of kp^GL is slightly smaller than its BEG value, 
however throughout the BGS-BEG evolution, kpXcL 
does not change substantially. 

The coherence length can be extracted from Eq. (jH) as 
= h/^2m\ a{T)\ lea ding to ^(T) = ^gMTT^/^ 
where £,gl ~ Tiyjc/ {2maa). Unlike the penetration depth, 
the coherence length ^gl changes substantially during 
the BGS to BEG evolution. In the BCS regime, kp(,GL ~ 
A/7C(3)/(127r2)(e^/Tc) is very large, and in terms of 
kpas becomes kp^^oL — ^exp[7r/(2fcF|as|)], where A — 
A/7C(3)/(127r2)(7re2-T/8). In the BEG regime, kp^GL = 



■\/7r/4 {\^\/ ep)^^^ is also very large, and in terms of kpag 
becomes kpS^Q^ 



'7r/4/-\/fcpa7. However, kp^cL passes 
through a minimun in the intermediate regime, where 

kp^GL « Oil). 

As discovered by Abrikosov [J], the parameter n = 
A(T)/^(T) is of fundamental importance in the charac- 
terization of the magnetic properties of charged super- 
fluids. When K < Kc = l/'s/S the charged superfluid is 
a perfect diamagnet (type-I), which does not allow the 
magnetic field to penetrate. When k > Kc, the charged 
superfluid allows the penetration of magnetic field in the 
superfluid state in the form of vortices (type-II). Since 
the temperature dependence of A(r) and £,{T) is exactly 
the same, the parameter 



A(r) A, 



GL 



mb 



BTTTqh^C^ 



(8) 



is independent of temperature. Notice that Ky^kpVq = 
s/ (167r), where s = kpb/{epc^) is a dimension- 
less parameter which is a function of l/{kpas) only. 
In Fig. [21 we show the evolution of func- 
tion of the scattering parameter l/{kpas). In the 
BGS limit, Ky/t^ = ^97r3/[14C(3)](rc/ei.) cor- 
responds to K^/kpTq = B exp[— 7r/(2A:i?|as|)], where 
B = ^97r3/ [14C(3)](8e''-V7r). While in the BEG limit, 
Ky^kpVq = \/2 {ep/\fi\)^^^ corresponds to n^kprq = 
\/2y/kpas. Notice the maximum of k in the vicinity of 
unitarity and = {l/[kpas] = 0.554). 

Unfortunately, in current experiments for neutral 
atoms in fictitious fields, only the product qH is con- 
trolled, instead of H alone [2|. So it is useful to think 
of ultra-cold superfiuids like ^Li or "'"K as having charge 
q ^ 0, but with qH finite. In this sense, these neutral 
superfiuids are extreme type-II with k — > cx3 throughout 
the BCS to BEG evolution 
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FIG. 2. Universal plot of the Ginzburg-Landau parameter k 
versus scattering parameter l/ikpas), where is the classical 
radius of a fermion with mass m and charge q. 

To obtain the phase diagram shown in Fig. [1] we set 
K = Kc in Eq. (j8|) and extract the fermion density ri as a 
function of l/{kpas), which leads to nr^ = s3/(15367r^). 
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In Fig. [U K is higher (lower) than Kc below (above) 
the critical line indicating a type-II (type-I) charged su- 
perfluid phase. For fixed density n, a phase transi- 
tion from type-I to type-II charged superfluid occurs, 
as the interaction parameter l/{kpas) increases. Elec- 
tron superfluids with lO'^^ cm~'^ < n < 10^^ cm~^ 
have 2.24 x lO'^^ < nr^ < 2.24 x lO'^s^ and the 
transition between type-I and typc-II occurs in the in- 
terval —8 < l/{kpas) < —4, while proton superfluids 
in nuclear matter, with w 3.60 x 10~** cm""^ and 
10^'' cm~^ < n < 10^® cm~^, have a type-I to type-II 
transition in the range —4 < l/{kpas) < 0. 

For type-I charged superfluids there is only the thermo- 
dynamic critical field Hc{T) determined by the condition 
HI{T)/{%tt) ^Fn- Fs, where F^ (F,) is the Helmhohz 
free energy for the normal (superfiuid) state. For a uni- 
form superfluid state the energy difference is Fn — Fg = 
aV(2/3) = a?/{2h) leading to H,{T) = i/c,o|e(T)| - 
|a(T)|-\/47r/6, where iJc,o = ao-\/47r/6 is the temperature 
independent prefactor. We define the dimensionless ther- 
modynamic critical field = Hcfi / Hkp , where Hkp ~ 
^qkp. Notice that -ffc.o = ^'^c,o/(27rei?), where Wc.o = 
\<l\Hcfl/{mco) is the cyclotron frequency at i?c,o- Using 
the asymptotic expressions for gq and b, we obtain i?c,o = 
yfc^V4/[77rC(3)](rc/eF) in the BCS regime, which 
can be rewritten as i?c,o = y^F^C exp [— 7r/(2fci?|as|)] , 
with C = ^4/[77rC(3)](8e'^-V7r). While we obtain 
Hc,o = V'^^(l/7r2)(ej./|/i|)i/4 in the BEC regime, 
which can be rewritten as Hc^o = \JkFrq{\/T:'^)\/kpas. 
The field i/c.o reaches a maximum near unitary and 
/i = 0, thus indicating that typc-I superfluids are most 
robust to the penetration of magnetic fields in that same 
region. 

For type-II superfluids there are two critical fields. 
The first is called Hc^{T) and separates the perfect- 
diamagnet Meissner phase from the non-uniform phase 
exhibiting vortices. The second is called i?c2 {T) and sep- 
arates the non-uniform phase exhibiting vortices from 
the normal state. Since ^Li and "'"K are extreme 
type-II superfluids with k — oo, then Hc^iT) — > 0, 
and thus we concentrate on the results for Hc2{T). 
The calculation of (T) is performed by lineariz- 
ing Eq. dH) -h^ {\/ -i2TrA/^qfi: + 2ma{T)i; = 0. 
Using the Landau gauge A = Hxy, the momen- 
tum components ky and k^ are good quantum num- 
bers and the solution for 'ip becomes ipn,ky,k-{x,y, z) = 
g(2fcBy+*'=2^)u„(2;), which substituted in the previous 
equation leads to the one-dimensional Schrddinger 
equation [—fi^ / {2'm)(fi / dx'^ + mujg^x — xo)^/2] Un{x) = 
e„Un{x), where xq = ^qky/{2nH) is the equilibrium po- 
sition of the harmonic potential, ujg ~ 2\q\H/{mca) is 
the harmonic potential frequency, and e„ — \a{T)\ — 
fi^k1/(2m) = huis{n + 1/2) is the eigenvalue. The 
highest magnetic field at which superconductivity nu- 



cleates occurs for n = and fcz = leading to 
the condition |a(T)| — TujJs/2. Isolating the mag- 
netic field from the harmonic potential frequency leads 
to H^.JT) = [^q/{2TT)]2m\a{T)\/h^, which can be fi- 
nally expressed in terms of the coherence length ^(T) 
as H^^iT) = $,/ [27rC2(r)] . Substituting £,{T) = 
eGLk(T)r'/', we write H,,{T) = H,,^oHT)l where 
^c2,o = ^q/i'^'^^GL)- Using again the reference field 
= ^qkp, we obtain the dimensionless upper critical 
field Hc2,o = Hc2fi/F[kp = l/(27rfc|,^Q^). This expression 
is equivalent to the ratio ?icjc2,o/(27reF), where cjc2,o = 
qHc2fl/i'mcQ) is the cyclotron frequency at i?c2.o- In the 
BCS limit, Hc2fi = [67r/7C(3)] {Tc/epf , which in terms 
of the scattering parameter l/{kpas) becomes Hc2,o = 
L>exp[-7r/(fcF|as|)] with D = 256e2T-V[77rC(3)]. In the 
BEC limit, Hc2fi ~ (2/^^) v^f/H which in terms 
of the scattering parameter \l{kpas) becomes i/c2,o = 
{2/iT'^)kpas. Since kpS^oL reaches a minimum in the re- 
gion near unitarity and /i = 0, it is clear that i?c2.o bas 
a maximum there, where type-II superfiuids are most ro- 
bust to the presence of real or fictitious magnetic fields. 

Before concluding, we note that the quantum regime 
htJc > 2ttT (where Landau level quantization is impor- 
tant) can be reached experimentally for ^Li and ''"K while 
preserving superfluidity. 

In conclusion, we have analyzed the effects of real or 
fictitious magnetic fields during the BCS to BEC evo- 
lution of s-wave superfluids with direct application to 
ultra-cold fermionic atoms. We have shown that a transi- 
tion from type-I to typc-II charged superfluidity occurs as 
the Ginzburg-Landau paramater crosses its critical value 
Kc ~ 1/V2 in the density versus interaction phase di- 
agram of fermions of charge q and mass m. We have 
shown that ®Li and in fictitious magnetic fields are 
extreme type-II superfluids. Finally, we have indicated 
that the critical magnetic flelds (real or fictitious) depend 
strongly on the scattering parameter l/kpUs and reach a 
maximum in a region near unitarity, where superfluidity 
is more robust to their penetration. 
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